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Abstract 

It is shown that any two Hamiltonians H(t) and H'(t) of ./V dimen- 
sional quantum systems can be related by means of time-dependent 
canonical transformations (CT). The dynamical symmetry group of 
system with Hamiltonian H (t) coincides with the invariance group of 
H(t). Quadratic Hamiltonians can be diagonalized by means of linear 
time-dependent CT. The diagonalization can be explicitly carried out 
in the case of stationary and some nonstationary quadratic H. Linear 
CT can diagonalize the uncertainty matrix cr(p) for canonical variables 
Pk, qj in any state p, i.e., a(p) is symplectically congruent to a diag- 
onal uncertainty matrix. For multimode squeezed canonical coherent 
states (CCS) and squeezed Fock states with equal photon numbers 
in each mode a is symplectic itself. It is proved that the multimode 
Robertson uncertainty relation is minimized only in squeezed CCS. 



1. Introduction 

The method of canonical transformations (CT) proved to be a fruit- 
ful approach in treating quantum systems. It is most efficient for systems 
which are described by Hamiltonians H which are quadratic in coordinates 
and momenta, or equivalently in boson creation and annihilation operators 
(quadratic systems) [1-6]. The main advantage of the method of CT consists 
in reducing the Hamiltonian H of the treated system S to a Hamiltonian H' 
of some simple system S' with known solutions. 
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Quadratic Hamiltonians H quad model many quantum (and classical) sys- 
tems: from free particle and free electromagnetic field to the waves in non- 
linear media, molecular dynamics and gravitational waveguide [7-10]. A con- 
siderable attention to nonstationary quadratic classical and/or quantum sys- 
tems is paid in the literature for a long period of time (see, for example, 
references in the review papers [9-12] and recent articles [13-15]. In the last 
decade they are intensively used in quantum optics to describe the nonclas- 
sical properties of light in squeezed states (SS) [7,16]. In these (and many 
other) fields the nonstationary H quad (t) play important role. 

Time dependent CT for quadratic systems <S qU ad (classical and/or quan- 
tum) have been considered in many papers [5,14,17-21]. The method devel- 
oped in [5] and intensively used in many papers (see references in [10]) con- 
sists of construction invariant coordinates and momenta qj(t),p^.(t), j,k = 
1, 2, . . . , N, [Qj(t),Pk(t)] = ihSkj (or invariant boson operators A^t), A ^ (t)) 
by means of time dependent CT. The generating operator for CT p^, qj — > 
Qj(t) 1S J us t the evolution operator of iS> qU ad- It is worth noting here that 
the new Hamiltonian H f , when expressed in terms of the invariant coordi- 
nates and momenta qj,p® is not diagonal, i.e. such CT are not diagonalizing. 

A CT is called diagonalizing if the new Hamiltonian H' in terms of the 
coordinates and momenta is diagonal quadratic form with constant (and 
positive) coefficients, i.e. H' is a Hamiltonian for a system of uncoupled 
harmonic oscillators (HO) H ho . One has to distinguish between two different 
types of diagonalization of H: (a) when H' is diagonal in terms of the new 
coordinates and momenta q'j,p' k , and (b) when H' is diagonal in terms of 
the old variables qj,Pk- If H is brought to H' by means of time-independent 
CT, then the spectrums of H and H' are the same since H and H' (and 
the corresponding states) are unitary equivalent. If the relating CT is time- 
dependent then H and H' are no more unitary equivalent (but the advantage 
of the unitary equivalence of the states remains). So we need a more general 
term of canonically equivalent Hamiltonians: H and H' should be called 
canonically equivalent if there is a CT, which transforms H into H' . 

Not any quadratic boson Hamiltonian can be diagonalized by means of 
linear time-independent CT [22-24]. Time dependent CT are much more 
powerful as we shall see below. Diagonalization of quadratic classical Hamil- 
tonian by means of time-dependent CT is studied by many authors [17-20]. 
In [18] the general N dimensional classical H quad (t) has been diagonalized us- 
ing linear time-dependent CT. In quantum case the diagonalization problem 
(the second type diagonalization) of the one dimensional HO with varying 
frequency uj(t) has been considered recently by Seleznyova [14]. 

The aim of the present paper is to consider the canonical equivalence of N 
dimensional quantum systems, in greater detail treating the case of systems 
with quadratic Hamiltonians. As applications we consider some general prop- 
erties of the dispersion matrix (called also fluctuation or uncertainty matrix) 
and multimode SS. 

In section II we show that any two N dimensional quantum Hamiltonians 
H(t) and H' can be canonically related via time dependent unitary operator 
U(t). In the case of quadratic Hamiltonians, H = H quad (t), the operator 
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U (t) is an exponent of a quadratic form of coordinates qj and momenta pk 
(that is, a general element the methaplectic group Mp(N, R)). In particular, 
such operators can diagonalize any H qVLa ^(t). We perform (in section III) 
the diagonalization of N dimensional H quad (t), expressing the parameters 
of the corresponding linear CT in terms of solutions of linear first order 
differential equations. For N — 1 these equations are reduced to the equation 
z + Q 2 (t)z = of classical oscillator with varying frequency. In section IV 
the diagonalization of the uncertainty matrix o for qj in any state by 
means of linear CT is established. Thus, a is symplectically congruent to 
a diagonal dispersion matrix. For squeezed canonical coherent states (CCS) 
[26] and for squeezed Fock states with equal photon numbers in every mode 
a is symplectic itself. 

2. Canonical equivalence of quantum systems 

The main aim of the method of CT is to reduce the Hamiltonian H of 
the treated system S to a Hamiltonian H' of some simple system S' with 
known solutions. CT in quantum theory are generated by unitary operators 
U, which we call the generator of CT. If CT is time-independent then H 
and H' are unitary equivalent and their spectrums are the same. However 
not any pair H and H' can be related by means of time-independent CT. In 
particular, not any quadratic H can be reduced to H' of harmonic oscillator 
by means of time- independent CT [22-24]. The time- dependent CT are much 
more powerful as we shall see below. 

Let \^{t)) be a solution of the Schrodinger equation [ihd/dt — H]\^(t)) = 
0. Then for any unitary operator U(t) the transformed state l^'^))) 

mt)) = u(t)Mt))- (i) 

obeys the equation [ihd/dt — H']\^'(t)) = with the new Hamiltonian H', 

H' = U{t)HU ] {t) - iW{t)dU ] {t)/dt. (la) 

Conversely, if two Hamiltonians H and H' are related by means of an (uni- 
tary) operator U(t) in accordance with eq. (la) then any solution \^(t)) of 
the system S is mapped into a solution \^'(t)) of the system S'. However, 
not any two given solutions \^(t)) and \^'(t)) of the two systems could be 
mapped into each other by means of U(t) since U(t) in general cannot act 
transitively in the Hilbert space. A more compact form of relation (la) is 
D'(t) = U(t)D(t)W(t), where D(t) = ihd/dt - H(t) (U(t) is interweaving 
operator for D(t) and D'(t)). From the requirement for the mean values of 
the "old" operator A and the "new" one A', 

(*(0|A|*(0> = (*'(*) I ^ I *'(*)>, (2) 
it follows that the operators A and A' are related as 

A! = U(t)AU\t). (2a) 
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Therefore the new canonical operators of the coordinates and momenta q' k 
and p' k , k — 1, 2, . . . , TV are related to the old ones as 

q' k = U(t)q k U^(t),p' k = U(t)p k U\t). (2b) 

Two quantum systems should be called canonically equivalent if there exists 
unitary operator U(t) which relates their Hamiltonian operators H and H' 
in accordance with eq. (la). U(t) could be called H-H' canonical equiva- 
lence operator (CEO). Let us note the main three advantages of establishing 
canonical equivalence of two systems (see also [14], where in fact canonical 
equivalence of one dimensional oscillators with constant and time-dependent 
frequencies was considered): 

(a) If we know solutions for one of the two canonically related systems we 
can obtain solutions for the other one by means of eq. (1). For a given system 
with Hamiltonian H the aim is to find simple H' with known solutions. As 
we show below this is always possible: any H in principle can be brought to 
as simple H' as that of a system of stationary oscillators or free particles. 

(b) If a time-dependent state \^>'(t)) of the system S' is an eigenstate of 
an operator A' then the U(t) partner \^{t)) = W(t)\^'(t)) in the system S 
is an eigenstate of the operator A = W(t)A'U(t). 

(c) If the operator A' is an integral of motion for S', that is A' commutes 
with the Schrodinger operator, 

dA'/Ot- (i/h)[A',H'] = 0, (3a) 

then the operator 

A = U ] (t)A'U(t) (36) 

is an integral of motion for the old system S. This property is very important 
since if we know one solution of a given S we can construct new solutions 
acting by the invariant operators on the known solution [5]. 
We shall prove now the 

Proposition 1. Any two N dimensional quantum Hamiltonians H and H' 
are canonically equivalent. The unitary operator U(t), which relates H and 
H' takes the form 



U(t) = Texp 



to 



H\t)dt 



UqT exp 



h 



f* H(t)dt] = S\t)U S\t), (4) 

J to 



where Uq is constant unitary operator and T and T stand for chronological 
and antichronological product. The solution (4) is unique for any initial 
condition U(0) = Uq. 

Proof. Let us perform two successive time-dependent CT by means of 
t/i = U S\t) and U 2 = S'(t), 



S*(t) =Texp H(t)dt 



(4a) 



Then from eq. (la) (taking into account dU{(t)/dt = (—i/K)HU\) we easily 
get Hi — for any Uq. The second transformation by means of U 2 = S'(t), 
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S'(t) = Texp ~j* H\t)dt 
then yields the required result: 



(46) 



H 2 = UHiW - ihU 2 dU\/dt = -ihS'dS'^/dt = H' . (4c) 

Now we see that the direct CT: H — > H' is performed by the operator (4) 
(note that the product of two unitary operators is also an unitary operator). 

The definition (la) of CEO U(t) can be safely multiplied by constant 
phase factors. We can prove that U(t), eq. (4), is uniquely determined by 
any initial condition U(0) = Uq. Indeed, suppose there is another unitary 
operator U(t) which also relates H and H' canonically and U(0) = U . Now 
we note that (it is easily derived from (la)) if U transforms H into H' then 
W transforms H' back into H and therefor the product V = ffiU keeps H 
invariant: 

H = VHV ] + l -[dV/dt]V ] and V(0) = 1. 

The latter can be rewritten in the form (3a) as dV/dt — (i/h)[V, H\ = 0, 
which means that V is an integral of motion for the system S. Any invariant 
operator for H has the form (note that S(t) is the evolution operator for S) 

V(t) = S(t)V(0)S*(t), (5) 

and since V(0) = 1 we have V(t) = S(t)S^(t) = 1. In a similar way one 
can get U(t)W(t) = 1 which proves that the two operators U(t) and U(t) 
coincide. This ends the proof of proposition 1. 

Thus, if initially two vectors ^'(O)) and |^(0)) are related by means of 
some Uo, 

|*'(0)> = £/o|*(0)>. (6) 

then at later time their evolutions \^'(t)) and \^(t)) governed by H' and H 
respectively are canonically related by means of U(t), eq. (4), with U(0) = 
Uq. If one takes |^'(0)) = |^(0)) then U(0) = 1 and one obtains unique 
U(t), eq. (4), which maps isomorphically the time evolved state of 
S on \^'(t)) of S' for all initial |^(0)). Note, eq. (6) doesn't imply that 
H(0) = H'(0): we have 

#'(0) = u q huI - ihu(o)[duydt] ]t=0 . 

Suppose now that H(t) and H'{t) are elements of a Lie algebra C Then 
S G G 3 S' , where G is the Lie group generated by L [27]. Thus, the CEO 
U(t) G G (for Uo — 1 and for Uq G G as well) and one can use the known 
properties of G [27] to represent U (t) in other factorized forms. 

The operator (4) converts canonically any iV dimensional H into any 
desired iV dimensional H' . In particular H can be converted into H' for a 
system of iV free particles or for a system of uncoupled harmonic oscillators 
(N mode free boson field). In the latter case if if is a quadratic form in terms 
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of N canonical operators q k and pj the operator (4) solves the diagonalization 
problem for quadratic Hamiltonians. 

We stress that one has to distinguish between two type of diagonaliza- 
tion: (a) H' is diagonal in terms of the new variables q' k p' k , and (b) H' is 
diagonal in terms of the old variables q k p k - In the first case the two systems 
S and S' respectively are treated in two different (q- and q'- ) coordinate 
representations (wave functions V(q,t) = (q\^(t)) and W(q',t) = (q'\^'(t))), 
while in the second case we can work in the same g-representation (wave 
functions V(q,t) = (q\V(t)) and W(q,t) = (q\W(t))). The second kind di- 
agonalization is achieved by means of operator U(t), eq. (4), with H' of the 
form of Hamiltonian H ho of N attractive oscillators, 
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N 



2 k 
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H' = -J2 —Pl + mkufal =H ho , (7) 



Let us briefly discuss the two CT, generated by Ui = UoS'(t) and U2(t) = 
S'(t) (see proposition 1). The first one brings H to zero, therefore the new 
states are time-independent. This is because S\t) is an evolution oper- 
ator for the S backward in time. After this CT the new canonical variables 

q^ = Ui(t)q k Ul(t) and pf = U^p.U^t) 
obey the equations (dUjdt = iU^H, dU\/dt = -iHU\) 

i.e., q k \ p^p are Heisenberg operators for the system S. The operator S'(t) 
is the evolution operator forward in time for the target system S'. 

It is worth noting here the CT generated by the evolution operator S(t) 
of S. This CT converts H(t) into H"(t) = S(t)H(t)S^(t) + H{t). If H is 
time-independent then H" = 2H. From (2) we derive that the new canonical 
variables after CT (9a), 

ql = S(t)q k S*(t) = q° k , p" k = S(t)p k S\t) = p° k , (9a) 

when expressed in terms of the old ones q k , pj, are integrals of motion of S, 
satisfying the eqs. (6a,b). Then the operators 

q°' k = U{t)qlU\t), p\ = U{t)plU\t) (9b) 

where U(t) is given by eq. (4), are integrals of motion for the system H' 
(recall the above property (c) of the canonical equivalence). Such integrals 
of motions for quadratic systems H(t) have been constructed in [5] and in- 
tensively used later [10]. 

If the target Hamiltonian H' is a sum of stationary attractive harmonic 
oscillators H-^o, eq. (7), in terms (for example) of the intermediate variables 
q^\ p k l \ Then the second CT (generated by ^(i) = exp(— (i/h)H\ 10 t) = 
Sho(t))- (q^ k \ P^) — ► WkiPk) yields the explicit relations between q k \ p^ 
and the final canonical variables q',p', 
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q k = q ( k ] cos(u k t) H p[ 1] sin(w fc i), (10a) 

P k = -rrikUJkq^ sin(uj k t) + p[ l) cos(uj k t). (106) 

This is an orthogonal symplectic transformation which for every k is a rota- 
tion in the plane q k \ (m k uj k )~ l p k lS> ). It is evident that the CT (9a) is not a 
diagonalizing CT, unless H itself is of the form H ho . 

It is of practical interest for a given Hamiltonian H to specify the class 
of CT (possibly time- dependent) which keep the H invariant (and thus they 
keep the Schrodinger equation invariant), i.e., 

H' = U{t)HU\t) - ihU{t)*^p- = H. (11) 

For time-independent U eq. (11) reduces to HU = UH. From (4) and (11) 
we have U(t) = S(t)UoS^(t), where Uq is arbitrary unitary operator. Thus, 
H-H CEO are unitary invariant operators. In the first paper of refs. [5] 
the dynamical symmetry group of a system S has been defined as a group 
of unitary operators which commute with ihd/dt — H and act irreducibly in 
the Hilbert space. Now we see that this symmetry group is highly nonunique 
since Uq is arbitrary and one can take it from irreducible representations of 
many groups. Then the set of invariants S(t)UoS^(t) realize an equivalent 
representation. For example, by means of the invariants q k and p° k one can 
construct an irreducible representation of (the Lie algebra of) the Heisenberg- 
Weyl group H W (N) and the quasi unitary group SU(N, 1) as well [5,27]. This 
means that the groups H W (N) and SU(N, 1) can be considered on equal as 
dynamical symmetry groups of any N dimensional system. 

In the next section we consider the above described canonical equivalence 
approach in greater detail for quadratic quantum systems, for which some 
explicit solutions can be obtained. 

3. Canonical transformations of quadratic systems 
and diagonalization. 

We consider the general N dimensional nonstationary quantum system 
with Hmiltonian H(t) which is a homogeneous quadratic form of coordinates 
and momenta (in this section H quaid = H), 



H(t) = A jk {t)pjp k + B jk (t)pjq k + B jk (t)qjPk + C jk {t)qjq k) (12a) 

where the coefficients Aj k (t) = A k j(t), Bj k (t), Bj k (t) and Cj k (t) = C k j(t) 
are arbitrary real functions of time. By no lack of generality one can put 
Bjk(t) = B k j(t). In (12a) the summation over the repeated indices is adopted. 
We can introduce iV component vectors 

Q= (9i,92, •••,9iv), P= (Pi,P2,-..,Piv) 



7 



and N x N real matrices A(t), B(t), C(t) (A(t) and C(t) are symmetric) and 
rewrite the Hamiltonian (12a) in a more compact form 



H(t) = pA(t)p + pB(t)q + qB(t) T p + qC(t)q, 

where B T is the transposed of B. To shorthand the notations it is con- 
venient to introduce the 2iV-vector Q = (p,q) and 2N x 2N matrix 7i (the 
grand matrix) and rewrite the Hamiltonian (12a) as (//, v — 1, 2, . . . , 2iV) 

H(t) = H^Q^Qv = QH(t)Q, (126) 

where U = ^ ^ T ^ . 

We note that nonhomogeneous quadratic Hamiltonians (i.e., Hamiltoni- 
ans (12a,b) with linear terms added) can be easily reduced to the form (12a,b) 
by means of simple time dependent displacement transformations. 

Let H' be an other quadratic Hamiltonian 



H'(t) = QH'(t)Q,H'=[ £r % . (13) 




Then the unitary operator U(t), eq. (4), which relates canonically Hamilto- 
nians (12) and (13), is an exponent of a quadratic in q and p form (we take 
U eMp(N,R)), 



U(t) = S'(t)U S*{t) = exp 



(14) 



where H(t) is a new grand matrix of the form (12b) and (13). H(t) can 
be expressed in terms of the Hamiltonian matrices 7~C{t) and 7i.'{t) using the 
Bakcr-Campbell-Hausdorff formula. In this case the operator (14) generates 
linear transformation of coordinates and momenta (we write it in iV x iV and 
2N x 2N matrix forms), 

Q'^m or (£) = (£ £)(f). («) 

where \ pp , \ pq , \ qp and \ qq are N x N matrices, 

A=( Xpp Xpq ) . (15a) 

From eqs. (la), (12b), (13) and (15) we obtain the following relation 
between the symmetric matrices 7i, 7i' and 7i (14) and the symplectic matrix 
A, 

i-H(t) = -n'(t) + A T H(t)A. (16) 

We see that for a given H(t) and 7i(t) this is a simple linear equation for 
H'(t). However for a given Hamiltonian matrices 7i and H' this is highly 
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nonlinear equation for 7i(t) since the matrix A(t) is to be expressed in terms 
of 7i(t) again: AQ = U(t)QW (t). Nevertheless for any given (differentiable 
with respect to t) matrices T~C(t) and 7~C'{t)) and for a given initial condition 
Tio the above equation has unique solution for 7i(t), since the expression of 
A in terms of Ti, is also differentiable and Peano theorem could be applied 
[37]. 

In this scheme A(t) is naturally represented as a product of two other 
2N x 2N matrices A^ and A^ 2 ^ of the form (15) corresponding to the two 
successive CT generated by U\(t) and f/ 2 (^) : 

A = A^A^; Q (1) = A«Q, Q' = A^Q^. (17) 

The matrices A^ 1 ) and A^ 2 ^ are seen to be solutions of the first order linear 
equations, 

i. A (i) = A (i) F (i)( t)) i. A (2) = f^A^, (18) 

(JjL (JjL 

where 

F«(i) = -2JH(t), FW(t) = 2JH'(t), J = ( °_i J ) > ( 18fl ) 
or in terms of the N x N matrices A, £>, C and A, £>', C, 

F<1, W = 2 (^ T 6 C )' ^W = -2(/ T /)■ (186) 

If H' is diagonal as for the oscillator system (7) then the second eq. (18) is 
easily solved: A^ 2 (t) = exp(2JHho^)Ao 2 ' ) . To perform the diagonalization of a 
quadratic H one has also to solve the first matrix equation in (18) and obtain 
A^(t), which in principle is always possible. In the case of stationary initial 
H the T exponent becomes ordinary one, so the explicit solution is given by 
the matrix exponent Aq 1 '* exp(— 2JH.t). The new variables are stationary H 
the diagonalizing CT is (highly nonunique) 

Q' = exp(2JH ho t)A ( o ) A^ ) exp{-2JH)Q. (18c) 

One can put Aq = 1, which corresponds to Uq in eq. (4). Having obtained 
explicitly A(t) = A^A^ one can next try to solve eq. (16) and obtain the 
generating operator U(t) in the form of the quadratic exponent (14). 

Note, the resulting H' is diagonal in the variables, which we choose for 
ifho- Let those variables be Then the final variables p' k , q' k obey 

eqs. (10a,b). Inverting the transformations (10a,b) we obtain H' diagonal in 
terms of the final variables as well: H' = Hh (ff, <f). In this way we perform 
explicitly the first kind diagonalization. If H' = iJho in terms of old variables 
Pk, Qk (second kind diagonalization), then H' is evidently not diagonal in 
terms of p' k , q' k . 

For some time-dependent H(t) explicit solutions of eqs. (18) can also 
be found. Thus, in the case of N — 1 following the scheme of refs. [4,5] 
one can express matrix elements of A^(t) (and therefore the CEO U(t) and 
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all quantum mechanical solutions of the system S) in terms of a complex 
function z(t) which obeys the equation of classical oscillator 

z + Q 2 (t)z = 0, 

where Q 2 (t) is simply determined by the parameters of the Hamiltonian (12), 

n 2 (t) = AAC + 2BA/A + A/(2A) - 3A 2 /(AA 2 ) - AB 2 - 2B. 

For HO with varying frequency uo(t) we have Vt 2 (t) = oo 2 (t). It is seen that 
an VL(t) corresponds to a class of quadratic H(t). For example constant Q 
corresponds to the stationary oscillator and to the oscillators with varying 
mass (damped oscillators) m(t) = m exp(— 2bt) and m(t) = m cos 2 bt con- 
sidered in [10,19]. Analytical solutions to the equation of z(t) are known for 
a variety of "frequencies" Q(t) [37]. In the case of an oscillator with varying 
frequency the diagonalizing CT generator U (t) has been expressed in terms 
oiz(t) in [14]. 

Let us briefly discuss the algebraic aspects of the above results. From the 
canonical commutation relations it follows that all 2N x 2N matrices A obey 
the relation (the symplectic conditions) 

AJA T = J, (19a) 

which for the N x N matrices \ qq , \ pp , \ qp and \ pq , defined in eq. (15)), read 

\q^pp ~ ^qp^pq = 1; ^qq^qp = ^qp^qqi ^pq^pp = ^PP^pq- (196) 

The set of matrices which obey the relation (19a) is defined as the symplectic 
matrix group Sp(N, R). It has N(2N + 1) real parameters. The rank of its 
Lie algebra is N. It is known that in classical mechanics the set of linear 
homogeneous CT forms a symplectic group Sp(N,R). In the quantum case 
the set of matrices A, which realize homogeneous linear transformations of 
the operators of coordinates and momenta close the same group. However the 
set of unitary operators U for which UqU^ and UpU^ are linear combinations 
of p and q contains one extra parameter, namely the phase factor. If one 
considers CT in greater detail as transformations of coordinates, momenta 
and vectors in Hilbert space one has to count the phase factors as well and 
then we get the larger group Sp(N,R) x (7(1) e Mp(N,R). If we consider 
transformations of coordinates, momenta and states we have to factorize over 
17(1): Mp(N, R)/U{1) = Mp(N, R). The resulting group Mp(N, R) is called 
methaplectic group. It is double covering of Sp(N,R). The Lie algebras 
of Mp(N,R) and Sp(N,R) are isomorphic [25]. They are of dimensions 
iV(2/A+l) and this is the number of independent matrix elements of matrix H 
in (14). The generators U (t) of linear CT (15) can be considered as operators 
of the unitary (but not faithful) representation U (A) of the symplectic group 
Sp(N,R). One can use the group representation technique [27] to represent 
U(t) G Sp(N, R) in several factorized forms. In the case of one dimensional 
nonstationary HO the CEO U(t) G SU (1,1) and its factorized forms have 
been considered in [14]. 
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— * 

If one considers Hamiltonians (12) with linear terms d(t)p + e(t)q added, 
then in the same way one would get that such inhomogeneous quadratic 
Hamiltonians can be diagonalized to the form (7) by means of the same 
U(t), eq. (4), this time U(t) being an element of the semidirect product 
group Mp(N, R)X) H W (N), where H W (N) is the N dimensional Heisenberg- 
Weyl group. 

4. Diagonalization of the dispersion matrix and 
multimode squeezed states 

The established possibility of converting (by means of time-dependent 
CT) any N dimensional Hamiltonian H to ifho of harmonic oscillators sug- 
gests to expect that the dispersion matrix a(Q,p) of canonical observables 
Q v , v = 1, 2, . . . 2N in any (generally mixed) quantum state p could be diag- 
onalized by means of some state dependent CT. It turns out that this really 
holds. 

Let us recall the notion of dispersion matrix a(Q,p) (called also fluctua- 
tion matrix, or uncertainty matrix). This is an 2N x 2N matrix constructed 
by means of the second momenta (the dispersions or the variances and co- 
variances) of coordinates q k and momenta pi (k, I = 1, 2, . . . , N; Q k = 
Pki Qn+u = Qk) 

<v = (y Q »Q„ = 2 (QnQ* + QvQn) - (Q») (Qv) , (20) 

where (Q) is the mean value of Q: for pure states we write (Q) = 
(fy | Q | and for mixed state p the mean values are (Q) = Tr (Qp). The 
matrix a is symmetric by construction. If we do not need to specify pure or 
mixed is the state (or the concrete choice of variable Q) we will simply write 
a(Q) (oi*(p)). 

For pure states |^), which are unitary equivalent either to multimode 
CCS \a) (a = (cci, a 2 , ■ ■ ■ , aw), N — 1, 2, . . .), or to Fock states \n) (n = 
(n 1 ,n 2 , • . . ,n N )) 

|^r) = U \a), or = U \n) (21) 

the dispersion matrix cr(Q, ^) can always be diagonalized by some CT. To 
show this let us note that for any CT, generated by some unitary U, 

Q^Q' = UQU\ (Q = (p,q)), (22) 

we have (using eqs. (20), (2) and (2a)) the following relation for the disper- 
sion matrices, 

a(Q', p) = a(Q, p'), a(Q>, tt) = <j(Q, (23) 

where the new state is p' = UpU^ and the new pure state is l^') = U\^f). 
Taking U = Uq and \a) = (or \n) = \^')) in eq. (23b) we obtain 

<r(Q', *) = a(Q, a) or <r(Q',^) = a(Q,n). (23b) 
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— * — * 

The dispersion matrices a(Q, a) and a(Q, n) are diagonal. Therefore, a(Q', ^) 
is also diagonal. This ends the proof that for any state |\&) which obey eq. 
(21) the dispersion matrix a(UQU\ ^) is diagonal. Note that a(Q,a) has 
the special diagonal form: 

a u (Q, a) = o- 22 (Q, «) = ... = cr NN (Q, a), (24a) 
a N+ljN+1 (Q, a) = cr N+2 ,N+2(Q, «) = ...= o- 2 n,2n(Q, a), (246) 

det a(Q, alf) = ]J a kk (Q, a) ■ a N+k ,N+k(Q, a) = ( —) N . (24c) 
k 4 

If \^{t)) is a time evolved state for a system S and |^(0)) satisfies the first 
equality in (21), then the time dependent CT (2) generated by U(t), eq. (4), 
with H' = H\y Q {p, q), and Uq taken as in eq. (21), will preserve in time the 
diagonal form (24) of a(Q'(t), ^(t)) since the evolution of CCS \a) governed 

— # 

by if ho is stable (up to a phase factor): \a) — > \a(t)}. 

All group-related CS \^(g)) [26,27] with initial (reference) vector |\P ) = 
| a) or l^o) — \n) (and any group G 3 g) by construction are unitary equiva- 
lent to | a) (or to \n)) , \^{g)) = f/o (^7) ^o) ? and thus their dispersion matrix 
obeys eq. (21), i.e. a(Q,^f(g)) is diagonalized to the form (24) by means of 
the (nonlinear) CT Q -> Q' = U (g)QU^(g). 

From this construction it is however not clear whether one can diagonal- 
ize a(Q,ty) for any state since in general \^) is not unitary equivalent 
neither to CCS \a) nor to Fock state \n). Surprisingly, it turned out that 
o~(Q,p) for arbitrary state p can always be diagonalized by means of linear 
CT. The crucial point is the positive definiteness of the dispersion matrix 
which is established by the following proposition (some known properties of 
the dispersion matrix are reviewed in, for example, ref. [10]). 

Proposition 2. The dispersion matrix a(Q,p) is positive definite for any 
(pure or mixed) state p. 

Proof. Consider the mean values of the positive definite z family of oper- 
ators F*(z)F(z), 

2N 

F(z) = zAQ = ]T z v AQ V , AQ„ = Q V - (Q„) , (256) 

V 

z v being arbitrary complex numbers (z is an 2N component complex vector). 
In any state p one has 

Tr[pF\z)F(z)} ee (F\z)F{z)) > 0. (26) 

The mean value (F^(z)F(z)) is expressed in terms of the dispersion matrix 
(20) as 

(F\z)F{z)) = z*[a(Q,p) - -TiJ]z = f $f > 0, (26a) 

where J is given by eq. (18a) and z*Jz = J2^,u z* Jy, v z v . The 2N x 2N 
matrix $ is nonnegative definite in virture of eq. (26). It is known [10] that 
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the dispersion matrix for canonical observables a(Q,p) is also nonnegative 
definite and obeys the Robertson inequality [28] 

det a(Q,p)>(^) N . (27) 

Now we shall prove that a(Q,p) is positive definite. In this purpose we put 
z = x + iy, x and y being arbitrary real vectors and write 

-> % -> -> 

< z*[a(Q,p) - -hJ]z — xa(Q,p)x + ya(Q,p)y- hxJy. (28a) 

If we put y = in (28a) keeping x arbitrary we would have xa(Q,p)x > 0, 
i.e. the dispersion matrix is nonnegative definite. From [29] we know that 
if a matrix a is nonnegative defined and its determinant is strictly positive, 
then a is positive definite. Thus, the dispersion matrix is positive definite, 

— * 

o~(Q,p)x > 0. End of the proof. 

Several interesting properties of the dispersion matrix o~(Q, p) follow from 
its positive definiteness. The first to note is related to its diagonalization. 
A well established fact is that any symmetric and positive definite 2N x 2N 
matrix a can be diagonalized by means of symmetric transformation with a 
symplectic matrix A [22-24]: AaA T = diag{d 1? d 2 , . . . , d 2A r}. Then we have 
the proposition 3, 

Proposition 3. The dispersion matrix a(Q,p) can be diagonalized by 
means of linear CT of coordinates and momenta, 



Q^Qd = ^dQ- o-(Q,p) -> cr(Qd,p) = diag{di,d 2 , . . . , d 2 jv}- (29) 
To prove the diagonalization formula (29) we have (in view of the positive 

— * 

definiteness of a(Q,p) and the result of ref. [22]) only to note the following 
transformation law of the dispersion matrix under linear CT Q' = AQ, 

a(Q',p)=Aa(Q,p)A T , (30) 

which is easily derived from eqs. (29) and the definition (20) of a . Note 
that: (a) the diagonalizing CT A d is not unique [25]; (b) A d is state depen- 
dent: Ad = Arf(p). Proposition 3 means that any quantum state is unitary 
and methaplectically equivalent to a state with vanishing covariances of co- 
ordinates and momenta. The general methaplectic operator (14) is called 
also squeeze operator [7,8,11,30], its canonical one dimensional form being 
exp[(^*a 2 — za) 2 )/2]. It is more adequate to call it squeeze and correlation op- 
erator since, e.g., for pure imaginary z it generates covariance of p and q and 
dosn't squeeze, while for real z it generates squeezing and doesn't correlate. 

The second property of o~(Q,p) which we note is that its determinant 
doesn't depend on the linear CT as it is seen from eq. (30), where det A = 1, 

det a(Q', p) = det a(AQ, p) = det a(Q, p). (31) 
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The trace of a is generally not invariant, however the quantities Tr[(aJ) 2fc ], k = 
1,2, ... , are invariant. Recall that time evolution governed by quadratic 
Hamiltonians is equivalent to linear CT of p, q. 

The third property of a(Q,p) to note is referred to the minimization of 
the Robertson inequality (27). We formulate it as a proposition. 

Proposition 4- The equality in the multimode Robertson uncertainty 
relation, eq.(27), holds in pure states |^) only and iff |\&) is unitary equivalent 
to a multimode CCS \a), 

\V) = U(A)\a), (32) 

where U(A) is an operator of the methaplectic group Mp(N, R). 

Proof. Let us first recall that Mp(N, R) is a quantum analog of the 
symplectic group Sp(N, R) 3 A and is generated by the set of operators 
{<lk<lu\{<lkPi + Piqk),PkPi}- Thus, U(A) in eq. (32) is an exponent of a 
quadratic form of coordinates and momenta, 

17(A) = exp[~(^HQ)], (32a) 

where H is an 2N x 2N matrix of the form (12b), related to A through the 
formula of CT AQ = UQlP . 

For further proof we need to diagonalize the dispersion matrix a(Q, p) by 
means of linear CT (29) and take into account eqs. (23a). We obtain 

a(Q d , p) = A d a(Q, P)Aj = <j(Q, p'),p' = T(A d )pT^(A d ), (33a) 

a(Q d ,p)=a(Q,p')=diag{d 1 (p'),d 2 (p'),...,d 2N (p')}, (336) 

where the diagonal elements d u (p') are the squared second momenta (the 
squared dispersions) of q k and p\ in the new state p' = T (A d ) pT^ (A d ) , 

dk(p') = o-kk(p') = Tr[(p fe - (p k )fp'} = a PkPk (p'), (34a) 

d N +k{p ) = o- N+k {p') = Tr[(g fc - {q k )fp] = v qkqk (p'), (346) 

Now we can write the determinant as a product of diagonal element d u (p'), 
v = l,2,...,2N, 

det a(Q,p) = det a(Q,p') = [d^d^p')] . . . [d N (p')d 2N (p% (35) 
From Heisenberg relation we have for every factor in eq. (35) the inequality 

h 2 

d k {p')d N+k {p') >—, k = 1, . . . , N. (36) 

From eqs. (35), (36) and (27) we derive that the equality in Robertson 
relation (27) holds iff the equality in eq. (36) holds for all modes (for every 
k = 1,2, ... N). And we know that this is possible iff p' is a pure state, 
namely a multimode CCS \a): 

p> ee U{A d )pU\A d ) = U(A d )\a)(a\U\A d ). (37) 
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This ends the proof of proposition 4. 

The Proposition 4 can be reformulated as follows: The equality in the 
multimode Robertson uncertainty relation (27) holds for pure states \^) only 
and iff \^/) is eigenstate of a set of boson operators b k , k — 1, 2, . . . , N, which 
are linear combination of the original a k and a k , a k = (2hm k u k )~ 1 ^ 2 (m k uj k q k + 
ipk), where m k and u k are mass and frequency parameters. 

Indeed, CCS \a) is an eigenstate of all a k with eigenvalues a k (ajfc|cf) = 
ak\6t)), and the minimizing states l^mm) are unitary equivalent to the CCS, 
l^min) = U(Ad)\a). Then, using (37) and the BCH formula, we see that 

^min 

) is eigenstate of b k , 

b k = U{A d )a k U ] {A d ) = u H ai + v kl a\, (38) 

feifcl^min) = a fe |* min ), (39) 

with the same eigenvalues a k . In view of (38) these eigenstates can be denoted 
as | a; u, v) or \a; A). 

For quadratic Hamiltonians the time evolution operator U qua a(t) G Mp(N, R). 
Therefore the time evolution of |^ m i n ) for quadratic Hamiltonians is stable, 
i.e., U(t)\fy min ) = \^' min (t)) up to a phase factor. \^' min (t)} is eigenstate of 
linear integrals of motions b k (t) with constant eigenvalues a k . Overcomplete 
system of such eigenstates \a;A(t)) (denoted as \a,t)) for N dimensional 
quadratic systems has been constructed in ref. [5] and used later in many 
papers [10]. In ref. [26] the states |0;A) = U(A)\0) (where |0) is annihi- 
lated by all boson operators a k ) are studied as Sp(N, R) CS with maximal 
symmetry and in ref. [30] the states |a;A) have been further studied and 
called multimode squeezed states. States |a;A) which minimize Robertson 
inequality should be called Robertson intelligent states (IS). They are of the 
form of squeezed CS. 

Let us note that in Hilbert space of pure states the representation U (A) 
doesn't act transitively and all states fall in Mp(N, R) orbits, which either 
coincide or have no common vectors. For example the two orbits U (A) \a) and 
U(A)\a') with a! ^ a are different. Now it is clear that due to the invariance 
of deta and Tr[(aJ) 2k ] under the linear CT, all states can be separated into 
larger classes with constant deter and Tr[(aJ) 2fc ]. For example the above two 
orbits U(A)\a) and U(A)\a') fall into one class with detcr(Q,^) = (h/A) N . 

A further property of the dispersion matrix a(Q,p) (the fourth one) we 
want to note here is referred to its symplectic character for a certain class 
of states. The matrix a(Q,p) is symplectic itself if it is symplectic for some 
canonical variables Q' . This follows from the transformation formula (30). 
When a(Q,p) is symplectic (more precisely it is a(Q, p) (det a(Q, p))~ l ^ 2N 
which is symplectic) it obeys the defining relation (19a), 

[a(Q,p)Ja T (Q,p)] [det a(p)]^ N = J, (40) 

since always det a (Q,p) > 1. In terms of the four N x N blocks o" pp (p), 
^(p), (T pq {p) and a qp (p) = cr^ q (p), 

V <W) MP) / V ^ ^4 J 
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the symplectic conditions (40) are rewritten as (I is N x N unit matrix) 

WM*M ~ VIM [^t v{p)r 1/N = /, (42a) 

0"4<T 2 — 0"Jo"4 = 0, <7i<73 — o\(7\ = 0. (426) 

Condition (42a) means that (if a(Q,p) is symplectic) the N x N matrix 
& qq (p)cr pp (p) — <Jp q (p) is a multiple of unity, the common factor being greater 
or equal to (h 2 /A) N (in view of relation (27)). 

We shall prove that relations (42a,b) hold in some of the N mode squeezed 
number states U(A)\n) (with det a > (h 2 /A) N if n 7^ 0) and in all squeezed 
CS (Robertson IS) \a; A) (with det a = (h 2 /A) N ). In the one dimensional case 
the states \a; A) = \a; u, v) are the Yuen two photon CS [32] and the relation 
(42a) is the equality in the Robertson-Schrodinger uncertainty relation: for 
N — 1 the quantities cr pp (p) and <r qq (p) are the squared variances of p and q 
and (Jpq(p) is their covariance. 

Before proceeding with the proof of (40) for the squeezed CS and Fock 
states let us note that one can use symplectic conditions (42a) to define 
a larger class of generalized multimode IS (or generalized multimode SS), 
which contains canonical IS \a; A) as a special case, corresponding to det a = 
(h/4) N . This can be called the class of SS with symplectic uncertainty matrix. 

If one consider the transformation properties of the N x N matrix com- 
binations occurring in (42a,b), o^O\ — o"^a 3 , 0\<J± — <7 2 a"J, <7 4 o" 2 ~ erf 1 04 
and <7iO"3 — crfVi, one could see that they are transforming linearly among 
themselves and the relations (42a,b) are invariant under linear CT. Then it 
is most simple to check the relations (42a,b) in the diagonal representation 
(33b) of the dispersion matrix, where they are reduced to one relation (since 
a 2 (Q,p') = = a 3 (Q,p')). 

a qq (p')a pp (p') / [det aip)} 1 ^ = I. (43a) 

In the diagonal representation det a{Q,p') = det[a qq (p')a pp (p')]. Eq. (43a) 
requires the N x N matrix a qq (p')a pp (p') to be a multiple of unity, i.e. 

a k (p')a N (p') = f(p') for all k = l,2,...,N. (436) 

If all products (Jkk(p')&N+k(p') are equal then the common factor f(p') must 
be invariant under linear CT and equal to [det ^"(p)] 1 ^ . 

— # 

<t(Q) is diagonal in CCS \a) and in Fock states \n), and eq. (43b) is 
satisfied in all CCS and in some Fock states. For Fock states the matrix 
a(Q,n) is diagonal with uncertainty products 

h 2 

<Jkk(Q, n)a N+k (Q, ft) = - (1 + 2n k ) 2 . (44) 

Thereby if n k are equal then the matrix a qq (Q,n)a pp (Q,n) is a multiple of 
unity with factor 

f(n) = (h 2 /A)(l + 2n k ) 2 = [det a(Q, n)} 1 ^ . (44a) 
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Thus, the two families of states U(A)\a) and U(A)\n) with nk = n obey eqs. 
(42a,b) and are states with symplectic dispersion matrix. Quadratic Hamil- 
tonians (and only they) will preserve this symplectic property in time evolu- 
tion. One can easily show that eq. (43b) can be also satisfied in states |^'), 
\^') = IIa: Ivlk, where [</?% are with vanishing covariance: & Pk q k (\<p') k) — 0. 

P.S. After the paper was completed (in August 1995) the author has learned 
about article [38], where the diagonalization of the dispersion matrix is established 
by means of several (explicit) linear CT. Some improvements in text are made in 
this e-print submission. 
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